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Abstract 

Umbral dots (UDs) are small isolated brightenings observed in sunspot um¬ 
brae. They are convective phenomena existing inside umbrae. UDs are usually 
divided into central UDs (CUDs) and peripheral UDs (PUDs) according to their 
positions inside an umbra. Our purpose is to investigate UD properties and 
analyze their relationships, and further to find whether or not the properties 
depend on umbral magnetic field strengths. Thus, we selected high-resolution 
TiO images of four active regions (ARs) taken under the best seeing conditions 
with the New Vacuum Solar Telescope in the Fuxian Solar Observatory of the 
Yunnan Astronomical Observatory, China. The four ARs (NOAA 11598, 11801, 
12158, and 12178) include six sunspots. A total of 1220 CUDs and 603 PUDs were 
identified. Meanwhile, the radial component of the vector magnetic field of the 
sunspots taken with the Helioseismic and Magnetic Imager on-board the Solar 
Dynamics Observatory was used to analyze relationships between UD properties 
and umbral magnetic field strengths. We find that diameters and lifetimes of UDs 
exhibit an increasing trend with the brightness, but velocities do not. Moreover, 
diameters, intensities, lifetimes and velocities depend on the surrounding mag¬ 
netic field. A CUD diameter was found larger, the CUD brighter, its lifetime 
longer, and its motion slower in a weak umbral magnetic field environment than 
in a strong one. 

Keywords: Sunspots, Umbra; Sunspots, Magnetic Fields 


1. Introduction 

There are some small bright features called umbral dots (UDs) in a dark umbra, 
and the features can be found almost all over umbrae and pores. UDs only 
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cover 3 “ 10 % of the umbral area, but contribute 10 ~ 20 % of the bright¬ 
ness. Therefore, convective motio ns must e xist within umbrae because radiation 
cann ot explain the phenomena (Deinzer, 1965 : Sobotka. Bonet. and Vazauezl . 
Il993h . Therefore, UDs play a vital role in the energy balance of sunspots. The 
study of UDs is crucial to understand convective motions and interactions of 
the plasma with strong magnetic fields, and analyze the formation mechanism 
of sunspots. Usually, UDs are divided into two clas ses according to their origins 
( Grossmann-Doerth. Schmidt, and Schrot^ 1986h . UDs are called peripheral 
UDs (PUDs) near an umbra-penumbra boundary and central UDs (CUDs) in an 
umbral center. PUDs are generally brighter than CUDs, and move quickly toward 
the umbral center, while CUDs are relatively static. Two different models have 
been proposed to explain the formation mechanism of UDs: the clustered mag¬ 
netic flux tube and the monolithic flux tube model. The former considers that 
UDs represeii t hot field-free gas intruded into a cluster of magnetic flux tubes 
(|Parkeilll979ll . The latter suggests that the energy transport in an umbra is dom¬ 
inated by non-stationary narrow rising plumes of hot plasma with adjacent down¬ 
flows. UDs are formed by the narrow up-flow plumes t hat become almost field- 
free near the surface layer ( Schiissler and Vbgleil 2006h . The essential difference 
between the two models is along the boundary of the UDs there exist local down¬ 
flows in the monolithic flux tube model. Detailed studies of UD properties, such 
as morphologies, velocities, lifetimes, and intensities, and relationships between 


different pro perties are crucial to understand th e nature of the local convec- 


tive motions dSobotka. B 

randt. and Simonl. 1997allbl: Bharti. Joshi. and Jaaffrev 

2007: Riethmiiller et a,l\. 

20081: Watanabe et aZ.|. |2012|l. 


Our purpose is to analyze UD properties and relationships between differ¬ 
ent properties, and further to find whether or not the properties depend on 
umbral magnetic field strengths. So we selected high-resolution observations of 
four active regions (ARs) taken under the best seeing conditions with the New 
Vacuum Solar Telescope (NVST) in the Fuxian Solar Observatory of the Yunnan 
Astronomical Observatory, China. 

The organization of the paper is as follows. The observations and data re¬ 
duction are described in Section 2. Section 3 briefly describes the identification 
procedure of CUDs and PUDs. Section 4 illustrates UD properties and their 
relationships. Moreover, the relationships between UD properties and umbral 
magnetic field strengths are discussed. Finally, we give our conclusions in Section 


5. 


2. Observations and Data Reduction 

The NVST is a vacuum solar telescope with a 985-mm clear aperture whose pur¬ 
poses are to obtain high-resolution imaging and spectral observations, including 
measurements of the solar magnetic field. The telescope consists of one channel 
for the chromosphere and two channels for the photosphere. The band used for 
observing the chromosphere is Ha (656.3±0.025 nm). The bands for observing 
the photosphere are TiO (705.8±1 nm) and G-band (430.0±0.8 nm). The high- 
resolution data of the NVST are classified into two levels. The level 1 data are 
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Figure 1. A reconstructed and projection-corrected image of NOAA 12158 recorded on 13 
September 2014 at 03:00:00 UT with the NVST. The active region included three sunspots 
marked as A, B, and C. 


processed by frame selection (lucky i maging;) dTubbj, 2004 ). The level 1+ data 
are reconstructed by speckl e masking dLohma nn. Weisrelt. and Wirnitzer . 1983) 
or iterative shift and add (Izhou and Lil. Il998l). Technic al details of the NVST 


were described by Lin et al. ( 2014h arid Xu et ai ( 2014h . 


Because the TiO line is highly sensit ive to umbral temperature variations 
(|Berdvngina. Solanki. and Frutigerl . 1200311 , its observations are more appropriate 
to investigate UDs. Therefore, we used the level 1+ TiO observations, rather 
than the Ha and G-band ones. We selected high-resolution image sequences of 
six sunspots taken under the best seeing conditions since October 2012. The 
sunspots located in four ARs: NOAA 11598, 11801, 12158, and 12178. The 
observation parameters are listed in Table 1. The high-resolution data were 
obtained without adaptive optics. But the resolution of the reconstructed images 
can almost reach as high a s the diffraction limit of the NVST in good seeing 
conditions ( Liu et oiH . 2014ll . From the table we see that the pixel size of the 
sequences is different before and after 2014. This is because the NVST team 
changed their optical system on 19 May 2014. 

The im ages in each sequence were co-align ed by the subpixel registration 
algorithm ( Feng et 'ai\ . 2012 ; Yang et 'oH . 12015 ). From Table 1 we see that the 
sunspots located in NOAA 11598, 11801, and 12158 were away from the so¬ 
lar disk center, especially the sunspot of NOAA 11801. Its heliocentric an¬ 
gle, 9, was close to 39° (ie., cos(0)=O.78). So we used the SolarSoft routine 
wcs-ConvertJrom_coord to transform the images of NOAA 11598, 11801, and 
12158 to heliographic coordinates for correcting the projection effects (ISun . 
2 OI 3 II . Figure 1 shows a reconstructed and projection-corrected image of NOAA 


12158 recorded on 13 September 2014 at 03:00:00 UT, where three sunspots (A, 
B, and C) are present. 
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Table 1. Parameters of the observed sunspots | 


AR NOAA 

Date (d-m-y) 

Time interval (UT) 

Location 

Pixel size (”) 

Cadence(s) 

11598 

29-10-2012 

05:50:13-07:44:36 

S11W27 

0.041 

37 

11801 

01-08-2013 

03:38:35-04:34:55 

W31N24 

0.041 

20 

12158 

13-09-2014 

02:57:30-03:47:14 

N20W27 

0.052 

30 

12178 

03-10-2014 

04:35:00-05:33:03 

S01E05 

0.052 

40 



3. Identification and Tracking of UDs 


We followed the method discussed in iFeng et al. (2015) to identify and track 
CUDs and PUDs. The method mainly consists of three steps: first, the periphery 
(umbra-penumbra) and center boundaries in an umbra are detected based on 
the morphological reconstruction technique; second, the UDs are identified based 
on the phase congruency technique; finally, the identified UDs are tracked based 
on a 26-connected neighborhood technique. The phase congruency technique has 
been us ed to extra c t low-c ontrast solar features, like coronal loops and umbral 
flashes ( Fener et I2m4allill . 

Empirical distance thresholds have been used to divide UDs into CUDs and 


PUD s ( Riethmullere|^^ , 200^ Watanabe. Kitai. and Ichimoto , 2009t Hamedivafa 
2011 : Louis et all . 2012ll . Riethmiiller et all ( 200^ considered an UD as a PUD 
if the UD’s birth position is closer t han 400 km f rom a defined umbra boundary, 
otherwise it is defined as a CUD. Hamedivafal (2011) defined a n arrow width 
near a n umbral boundary where the UDs are considered as PUDs. iLouis et all 
( 2 OI 2 II considered UDs with a threshold 0.8” inward from an assigned umbra 
boundary as CUDs, and those located outward as PUDs. 

Because the dehnition of the periphery (umbra-penumbra) and center bound¬ 
aries is crucial to classify UDs, we briefly intr o duce the dehnition and the 
identihcation procedure proposed bv IPeng et all (201^. They dehned the pe¬ 
riphery and center boundaries according to an umbral prohle. We show the 
three-dimensional surface of a reconstructed and projection-corrected image ob¬ 
tained from NOAA 11801 in Figure 2. The X and Y axes indicate the image size 
and the Z axis its normalized intensity. The color bar on the right indicates the 
intensity range from 0 (blue) to 1 (red). As illustrated with a red dashed line, the 
prohle of the umbra can be approximated by a trapezoid whose two sides appear 
skewed, and the base is relatively hat. Here, the base of the trapezoid is dehned 
as the center of the umbra where the UDs are considered as CUDs. The two 
skewed sides of the trapezoid are dehned as the peripheral region where the UDs 
are considered as PUDs. For obtaining the periphery and center boundaries, an 
image was reconstructed by the morphological reconstruction technique, and two 
thresholds, 0.3 and 0.6 Rmax, were used. Rmax denotes the maximum intensity 
of the reconstructed image. 
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Figure 2. The 3D surface of a reconstructed and projection-corrected spot in AR 11801 whose 
original image was recorded on 1 August 2013 at 03:48:04 UT. The X and Y axes indicate the 
image size and the Z axis its normalized intensity. The color bar on the right indicates the 
intensity range from 0 (blue) to 1 (red). The red dashed line illustrates the umbral profile that 
can be approximated by a trapezoid. 


Identified UDs were divided into CUDs and PUDs according to the bound¬ 
aries. However, we only obtained an insignificant number of PUDs in sunspots 
A, B, and C of NOAA 12158. The PUDs are so few that we cannot obtain 
accurate statistical results. Therefore, we abandoned the PUDs in these three 
sunspots. As a result, we extracted CUDs from six sunspots and PUDs from 
three sunspots. Although the division method of CUDs and PUDs might mis- 
classify a few CUDs or PUDs, this has no effect on statistical results of UDs. Two 
identified results are shown in Figure 3. The periphery and center boundaries 
are marked with yellow and red curves, and the UDs are marked with white 
contours. 


4. Results and Discussion 

Firstly, the means and standard deviations of UD properties such as equivalent 
diameters, ratios of the maximum intensity to the average intensity of the cor¬ 
responding adjacent umbral background, lifetimes, and horizontal velocities are 
obtained with distribution functions. Secondly, the intensity-diameter, lifetime- 
diameter, lifetime-intensity, and velocity-intensity relationships are analyzed. 
Finally, the relationships between UD properties and umbral magnetic held 
strengths are discussed. 

4.1. Property Dehnition and Feature Extraction 


For each UD, its equivalent diameter (D^g) is calculated as -^dA/Tr, where 
A denotes the total number of pixels of an UD. The ratio of the maximum 
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Figure 3. Boundaries identified by the method llFeng et al.l l2Q15fl . The periphery (um¬ 
bra-penumbra) and center boundaries are marked with yellow and red curves, and the UDs 
are marked with white contours, (a) The result for a corrected sample image from NO A A 
11801 whose original image was taken on 1 August 2013 at 03:48:04 UT. (b) The result for 
the sunspot marked as B in Figure 1. Note that there is no yellow curve. This is because only 
an insignificant number of PUDs were identified, and they are not enough to obtain accurate 
statistical results. Thus, we abandoned the PUDs and the corresponding peripheral region. 


Table 2. The number of the identified UDs of each sunspot 


AR NOAA 

11598 

11801 


12158 


12178 

Spot 



A 

B 

C 


CUD 

313 

201 

90 

279 

75 

262 

PUD 

219 

162 




222 


intensity (lud) to the average intensity of the adjacent unibral background (Ibg) 
are determined from the identified areas. Some UD properties are defined using 
the tracking procedure; they are lifetime (Tud)-, birth-death distance {Lm): and 
horizontal velocity (Vud)- Tud is the sum of the cadence of all frames; Lm is the 
centroid distance from its birth to death frame; Vud is computed by dividing L^d 
by Tud- 

We rejected the UDs whose diameters were lower than 130 km (0.18”) and 
lifetimes were less than two minutes for accurate statistical results. In the track¬ 
ing procedure, if splitting or merging occurred, the UD was discarded. As a 
result, a total of 1220 CUDs and 603 PUDs were identified. Table 2 lists the 
number of CUDs and PUDs of each sunspot. Subsequently, we utilize them to 
analyze UD properties and their relationships, and the relationships between 
these properties and umbral magnetic field strengths. 

4.2. UD Properties 

We obtained the probability histograms of UD properties in each sunspot, and 
fitted them with distribution functions for obtaining property values of UDs. 
Figure 4 shows the histograms and their fit curves of the UDs located in AR 
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Table 3. All means and standard deviations of UD properties 


CUDs 


NOAA 

Spot 

Diameter (km) 

^ud/^bg 

Lifetime(min) 

Velocity(km s 

11598 


178±40 

1.05±0.02 

5.35 

0.37±0.19 

11801 


216±38 

1.07±0.03 

4.48 

0.30±0.16 

12158 

A 

235±41 

1.10±0.05 

5.67 

0.27±0.14 

12158 

B 

225±40 

1.08±0.04 

5.66 

0.38±0.20 

12158 

C 

234±45 

1.09±0.04 

7.26 

0.19±0.10 

12178 


210±38 

1.07±0.03 

4.59 

0.37±0.19 

PUDs 

11598 


195±41 

1.08±0.06 

6.25 

0.47±0.24 

11801 


226±41 

1.12±0.05 

8.12 

0.51±0.27 

12178 


226±46 

1.15±0.09 

7.95 

0.45±0.24 


NOAA 12178; other ARs are not shown due to similar histograms and fit curves. 
But all the fitted means and standard deviations are listed in Table 3. The red 
color indicates the histograms and the corresponding fit curves of the CUDs, 
and the blue color those of the PUDs in Figure 4. 

Figure 4a and b show the probability histograms of equivalent d iameters and 
inten s hy ratios. We fo ll owed the method used i n previous works (IHamedivafa 


2011 : Watanabe et al. . 2012t iFeng et l2015[) to fit the distribution of the 


diameters and the intensity ratios with a Gaussian function. The function is 
as follows. 

1 


y = 






( 1 ) 


Here, ^ and cr are the mean and the standard deviation, respectively. ^ was 
used as the property values of equivalent diameters and intensity ratios. The 
mean diameters of CUDs range from 178 to 235 km, and that of PUDs from 195 
to 226 km (see the third column in Table 3). The diameters of all the UDs range 
from 1 31 to 366 km. The results are in good aKreement with those of previous 


work s (IH amedivafa. 200^ Riethmuller et ai] . l2008 : Bharti. Beeck. and Schusslcrl .[ 


201fll : Kilcik et al . 2012ll . They concluded that the range of UD diameters is 
between 150 and 350 km. The mean intensity ratios of CUDs range from 1.05 
to 1.10 and that of PUDs from 1.08 to 1.15 (see the fourth column in Table 3). 
This demonstrates that the brightness of CUDs is 5 - 10 % higher than that of 
their adjacent backgrounds, while the brightness excess is 8 - 15 % for PUDs. 

As shown in Figure 4c, the pr obability histograms of UD lifetimes approx i- 
mate an exponential distribution. Rlethmiiller et al. (2008); Hamedivafal (2011); 
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Figure 4. Probability histograms and fit curves of the UDs located in the sunspot of NOAA 
12178. (a) Equivalent diameter, (b) intensity ratio, (c) lifetime, and (d) horizontal velocity. 
The red color indicates the histograms and fit curves of CUDs, and the blue color those of 
PUDs. 


Fenef et al\ 1)2015fl also found that UD lifetimes exhibit a similar distribution. 


So we used an exponential function to fit the histograms, which is described as 
follows. 


y = Xe 


— \x 


( 2 ) 


For an exponential distribution, both the mean and the standard deviation 
are equal to A“^, which was used as the property value of UD lifetimes. The 
mean lifetimes of UDs range from 4 to 8 minutes (see the fifth column in Table 
3). This result is in qualitative agreement with other works in the literature 


(ISobotka et al . lOOOl: Riethmiiller et oL . 20081: Hamedivafa . 2011 : Feng et al 


201 ,^ . 


The horizontal velocity histograms of UDs, regardless of CUDs and PUDs, 
show approximately Gaussian distribution with zero mean and the same variance 
in the X and Y axes of the helioprojective-Cartesian coordinate syste m, demon¬ 
strat ing that the horizontal velocity follows a Rayleigh distribution ( Siddiouil . 
19641) . Therefore, we used a Rayleigh function to fit the velocity histograms. The 
function is dehned as. 


X _ 

y = 

rr^ 


( 3 ) 
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where a is the scale parameter of the distribution. The mean and the standard 
deviation are equal to ayJ'K 12 and — 7r)/2, respectively. We used the two 

values to describe the velocity property of UDs, rather than the peak of the 
distribution. The mean values and standard deviations are listed in column 6 of 
Table 3. The velocity range of CUDs is between 0.19 and 0.38 km and that 
of PUDs between 0.45 and 0.51 km (see Table 3); the horizontal velocity 
of CUDs is slow e r tha n that of PUDs. IWatanabe. Kitai. and Ichimotol ( 200911 
and iKilcik et al\ ( 2012ll concluded that PUDs have a relatively fast horizontal 
velocity because they are located in regions with a strong horizontal component 
of the magnetic field and/or a strongly inclined field. Our findings confirm their 
conclusions. 


4.3. Correlation Analysis of UD Properties 


The relationships, such as intensity-diameter, lifetime-diameter, lifetime-intensity,! 
and velocity-intensity, for all CUDs (1220, obtained from six sunspots) and all 
PUDs (603, from three sunspots) are presented in the form of scatter plots in 
Figure 5. 

For obtaining the trend of the scatter points, the points were sorted into 
15 equally spaced bins along the X axis between the minimum and maximum 
values of data. The average value of the data within a bin is indicated with a red 
square symbol, and the green vertical solid line denotes its standard error bar. 
The locations of the red squares specify the right side of each bin on the X axis, 
thus the maximum value of the data along the X axis is the location of the last 
square. A weighted linear fit to the average values of the bins is shown with a red 
dashed line in each panel of Figure 5. The weight of each bin is determined by 
the nu mber of the data in the bin. The method o f binning and linear fits has been 
used b y lWatanabe. Kitai. and Ichimotol ( 2009ll : iKilcik et al. ( 2012ll : Louis et ^ 

(I 2 m 2 ll . 

In the left column of Figure 5, the intensity-diameter (Figure 5a), lifetime- 
diameter (Figure 5c), and lifetime-intensity (Figure 5e) exhibit an increasing 
trend. However, the velocities present a weak inverse relationship to the intensity 
ratios (Figure 5g). The results demonstrate that larger CUDs tend to be brighter, 
live longer and move slower. The trends of PUDs are similar with CUDs, which 
can be found from the panels in the right column of Figure 5. The weighted 
co rrelation coefficients of the bi n s are listed in Table 4. 

Bharti. Beeck. and Schhsslei ( 20101) concluded from realistic radiative MHD 


simulations that the area-lifetime and brightness-lifetime relationships are posi¬ 
tively correlated. They stated the reason is that stronger and more extended con¬ 
vective up-flows are mai ntained longer and create larg e r and brighter UDs. Simi¬ 
lar tre n ds were found by Tritschler and Schmidt ( 2002 ); Sobotka and Puschmann^ 
( 2009ll : Kilcik ( 2012t) . Kilcik et air ( 2012h studied the statistical properties 
of UDs using high-resolution observations recorded by the New Solar Telescope 
at the Big Bear Solar Observatory and three-dimensional MHD simulations of 
sunspots. They concluded that the UD velocities are inversely related to their 
lifetimes. Our finding as regards of the CUD and PUD trends also supports their 
conclusions. 
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(a) 1 Intensity-Diameter relationship of CUDs 

320 

■ (b) 1 Intensity-Diameter relationship of PUDs| 
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Figure 5. Scatter plots of all UDs obtained from the four ARs. The left column illustrates 
the CUD relationships and the right column the PUD ones, (a—b) Corresponds to the intensi¬ 
ty-diameter relationships, (c-d) to the lifetime-diameter, (e-f) to the lifetime-intensity, and 
(g—h) to the velocity-intensity ones. The red square symbols indicate the average values of 
bins, and the green vertical solid lines denote their standard error bar. The red dashed lines 
are the weighted linear fit of the bins. 
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Table 4. Weighted correlation coefficients of UD properties | 



Intensity-Diameter 

Lifetime-Diameter 

Lifetime-Intensity 

Velocity-Intensity 

CUD 

0.95 

0.90 

0.89 

-0.40 

PUD 

0.95 

0.92 

0.97 

-0.03 



4.4. Relationships Between UD Properties and Umbral Magnetic Field Strengths 


In order to analyze relationships between UD properties and umbral magnetic 
held strengths, we used the radial component of the vector magnetic held, Br 
for the six sunspots taken with the Helioseismic and Magnetic Im ager on¬ 
board the Solar Dynamics Observatory iSDO/HMI. ISchou et a,l. . 2012h . Because 
the HMI image pixel- s cale i s 0.5” and UD diameters are approximately 0.3 


(ISchiissler and Vogleil 20061 : iKilcik et all l2012t iFeng et aj.. . 201,^ Yan et at 


201 ,^ . we failed to obtain an accurate pixel-by-pixel magnetic held feature for 
each UD. So we used the average magnetic held strength of the peripheral and 
central regions of each umbra. To improve the sensitivity of the average held 
strength, the maps were selected during the observed time interval of each 
data set and then averaged. Figure 6 shows two Br maps that were from NOAA 
11801 and 12158. The maps were hrst extended to the scale of TiO images 
using a nearest-neighbor interpolation method and aligned to the corresponding 
corrected TiO image (i.e., the image shown in Figure 3). In this hgure, the 
periphery and center boundaries are superposed on the maps with yellow and 
red curves. The average magnetic held strength of the peripheral and central 
regions are listed in Table 5. 

Figure 7 shows the relationships between CUD properties (diameter, intensity, 
lifetime, and velocity) and umbral magnetic held strengths. In each panel of 
Figure 7, the six red square symbols indicate the property values of six CUD 
data sets obtained with six sunspots, and those corresponding green vertical 
solid lines denote their standard error bars; the red dashed line represents the 
best linear ht of the red squares. These correlation coefficients are -0.87, -0.91, 
-0.47, and 0.79 for the four hts in Figure 7a-d. 

From the htted lines in Figure 7 we hnd that diameters, intensities and life¬ 
times of CUDs have a decreasing trend with increasing magnetic held strength 
(see plots in Figure 7a-c), while velocities increase (see Figure 7d). The results 
demonstrate that CUDs are larger and brighter, and their lifetimes longer, how¬ 
ever, their motions slower in a weaker umbral magnetic held environment than 
a stronger one. To PUDs, the three means of the properties obtained with three 
sunspots are so few that we did not obtain a signihcant result. 


5. Conclusions 

We selected high-resolution TiO image sequences of four ARs taken under the 
best seeing conditions with the NVST to investigate UD properties and analyze 
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Figure 6. The radial component of the vector magnetic field Br superposed on the periphery 
(umbra-penumbra) and center boundaries with yellow and red curves: (a) Br map of NOAA 
11801 obtained with the SDO/HMI on 1 August 2013 at 03:48:00 UT. (b) Br map of NOAA 
12158 obtained on 13 September 2014 at 03:00:00 UT. 

Table 5. Mean magnetic field strengths (G) of each umbra in the central 
(surrounding CUDs) and peripheral (surrounding PUDs) regions. 


CUDs 

AR NOAA 

11598 

11801 

12158 


12178 

Spot 



A B 

C 


Mean Field Strength 

2334 

2009 

1566 2079 

1741 

2212 

PUDs 

Mean Field Strength 

1888 

1582 



1688 


their relationships. Subsequently, umbral magnetic field strength relationships to 
UD properties were analyzed using the radial component of the vector magnetic 
field obtained from SDO/HMI. 

We found that diameters and lifetimes of CUDs hold an increasing trend with 
the brightness, but velocities do not. The trends of PUDs are similar to those of 
CUDs. Moreover, UD properties depend on their corresponding magnetic field 
environment. A CUD diameter is larger, its brightness higher, its lifetime longer, 
and its velocity slower in a weak umbral magnetic field environment than in a 
strong one. 
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Figure 7. Relationships between CUD properties and umbral magnetic field strengths: (a) 
diameter-magnetic field strength plot, (b) intensity-magnetic field strength plot, (c) lifetime— 
magnetic field strength plot, and (d) velocity-magnetic field strength plot. The red square 
symbols indicate the means of CUD properties in each sunspot, and the corresponding green 
vertical solid lines denote their standard error bars; the red dashed lines represent the best 
linear fit to the property values. 
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